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We report a thermally activated metastability in a GaAs double quantum dot exhibiting real-time 
charge switching in diamond shaped regions of the charge stability diagram. Accidental charge traps 
and sensor back action are excluded as the origin of the switching. We present an extension of the 
canonical double dot theory based on an intrinsic, thermal electron exchange process through the 
reservoirs, giving excellent agreement with the experiment. The electron spin is randomized by the 
exchange process, thus facilitating fast, gate-controlled spin initialization. At the same time, this 
process sets an intrinsic upper limit to the spin relaxation time. 

PACS numbers: 


Spins in quantum dots [I] are promising candidates 
for the realization of qubits - the elementary units of a 
quantum computer. Great progress was made in recent 
years towards implementing quantum information pro¬ 
cessing schemes with electron spins in GaAs quantum 
dots Hi], which hold the potential for scaling to a large 
number of qubits (SHU- Stable qubits with long coher¬ 
ence times are of crucial importance to execute numerous 
coherent quantum gates. Spin echo and dynamical de¬ 
coupling techniques were successfully employed to isolate 
the electronic system from the slowly fluctuating nuclear 
spins of the GaAs host material 0Q2HIE enhancing the 
coherence time X \ from below 1 /is to much longer times 
exceeding 0.2 ms. A fundamental limit < 2 T\ is set by 
the spin relaxation time T\. In a magnetic field, spins re¬ 
lax through spin-phonon coupling mediated by the spin- 
orbit interaction (3 HMEE]. Since here the spin-orbit 
coupling is weak, very long T\ times result, exceeding 1 s 
at 1 T [19], leaving ample room for further improvements 
of the spin qubit coherence. 

In this Letter, we report the experimental observation 
of a thermal electron exchange process via the reservoirs 
of a quantum dot, setting an intrinsic upper bound to T\ 
which can be orders of magnitude lower than the funda¬ 
mental spin-phonon limit m- The resulting metastable 
charge states - appearing in the double dot (DD) in ab¬ 
sence of interdot tunneling - make the exchange pro¬ 
cess detectable with a charge sensor. Within a diamond 
shaped region, the DD switches its charge-state back and 
forth over time from an electron on the left dot to an 
electron on the right dot without direct interdot tunnel¬ 
ing. After excluding unintentional charge traps and sen¬ 
sor back action, we present an extension of the orthodox 
DD transport theory accounting very well for the obser¬ 
vations. The exchange process can be used for fast qubit 
initialization [7]. Finally, we outline ways to extend T\ 
up to the spin-phonon limit. 

The sample is fabricated from a GaAs heterostructure 


with a 2D electron gas 110 nm below the surface (density 
2.6 • 10 11 cm -2 and mobility 4 • 10 5 cm 2 /Vs). The de¬ 
vice layout, see Fig.[lja), is adopted from Ref.[20| Each 
dot adjacent to the DD (center) acts as a charge sen¬ 
sor [21] 22] . changing conductance g up to a factor of 
two upon adding one electron to the dot closer to the 
sensor. The high sensitivity is due to the steep edges 
of a Coulomb blockade peak where the sensor is biased 
}20] . Simultaneously, strong capacitive shifts of the sen¬ 
sor biasing point are seen when changing DD gate volt¬ 
ages. Thus, compensation with linear feedback on sensor 
plunger L2 is employed to maintain charge sensitivity. 

The experiment was done in a dilution refrigerator at 
base temperature T ~ 20 mK. Ag-epoxy microwave fil- 
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FIG. 1: (a) Scanning electron micrograph of a nominally iden¬ 
tical device, with contacts (yellow) and gates (white boxes). 
Gates pl-p3, wl, wr and n are used to form the DD in the 
center while L1-L3 control the left and R1-R3 the right sensor 
dot. (b) Numerical derivative Ag of the left sensor conduc¬ 
tance. (c) Larger CSD with gate configuration as in Fig. [5] 
Circles indicate vertices where the metastability was observed. 
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FIG. 2: Sensor A g (a) and standard deviation of g (b) from 
20 ms measurement windows with reduced tunnel rates com¬ 
pared to Fig. B b ) . Dashed lines indicate the borders of the 
metastable diamond while solid lines show the reservoir tran¬ 
sitions, obtained from (a). The white and yellow dots show 
gate voltage configurations used in Figs. [3] and [4] 


ters and thermalizers [23 are mounted at the mixing 
chamber for improved cooling [231 25] « giving an elec¬ 
tron temperature of T e ~ 60 mK from Coulomb blockade 
thermometry [26]. All data presented here was acquired 
with the left sensor dot, though the right sensor gives 
essentially the same results. The DD reservoirs are held 
at the same potential. The sensor is typically biased at 
60 /jlV DC and the resulting current is digitized in real¬ 
time with a measurement bandwidth of 10 kHz, limited 
by the signal to noise ratio, not technical bandwidth. 

The charge stability diagram (CSD) of the DD is shown 
in Fig.[l][c). Sharp lines indicate charge transitions, form¬ 
ing the usual DD honeycombs E3- Charge states are la¬ 
beled (TVi, N 2 ) denoting the absolute number of electrons 
on the (left, right) dot in the ground state. The first two 
triple points are shown in Fig.JlJb). Dark lines designate 
transitions involving a reservoir, while the bright line cor¬ 
responds to equal left and right dot energies (zero detun¬ 
ing), separating the (1,0) and (0,1) states. The absence 
of curvature indicates weak interdot tunneling. These 
data reflect standard DD behavior as expected m- 

The CSD drastically changes upon reduction of the in¬ 
terdot and reservoir tunnel rates, see Fig. [2j The zero de¬ 
tuning line is no longer present, transforming instead into 
a broad region where the DD enters a metastable charge 
state, repeatedly switching between the (1,0) and (0,1) 
configurations over time, see Fig.[3](e). The switching is 
also visible in the standard deviation of the sensor sig¬ 
nal, see Fig. Bb), and fits very well into a diamond drawn 
with lines following the slopes of the reservoir transitions. 
This suggests involvement of the reservoirs in the switch¬ 
ing process. The measurement bandwidth is too low here 
to resolve the reservoir transitions in the standard devi¬ 
ation, thus decreasing the size of the diamond compared 
to a straight continuation of the reservoir lines. 

To quantify the switching rates, we gather numerous 
events throughout the diamond and histogram the dwell 
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FIG. 3: (a),(b) Switching rates on a log color scale, (c) 

Switching frequency / on a linear color scale, (d) Probability 
Pr of finding an electron in the right dot. (e) Real-time trace 
and (f) Y r (T) with exponential fit (red) to the high-T data, 
taken at the center of the diamond [white dot in Fig.[2jb)] . 


times in either charge state, finding single exponential 
decays. The rate Tl for switching into the left dot is given 
by the total number of switches N into the left dot and 
the accumulated waiting time t R in the right dot, Yl := 
F(o, 1)->(1,0) = N/tR. Similarly, V R is obtained, and Y rr 
are shown in Fig.[3|a,b). These rates are largest along the 
upper (a) or lower (b) edges, where a dot level approaches 
its reservoir level, here incidentally reaching the sensor 
bandwidth. Away from these edges, the rates fall off 
exponentially by three orders of magnitude until a small 
background rate is reached which varies weakly with gate 
voltages and is well above the smallest detectable rate. 
Photon assisted tunneling, residual interdot tunneling or 
a higher order process, possibly involving phonons m, 
could give rise to such a background rate. 

The switching frequency / = (T^ 1 + T^ 1 ) -1 of the 
complete process, i.e. from left to right dot and back, 
peaks close to the triple points, see Fig.[3](c) , and is low 
elsewhere. The probability Pr of finding an electron in 
the right dot is Pr = Yl/(Yl + Fr), shown in Fig.[3|d), 
reproducing the standard DD CSD with equal proba¬ 
bilities for an electron on the left or right dot at the 
zero detuning line. In addition, the measured temper- 
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ature dependence of Tr is shown in Fig.[3jf) , recorded 
at the center of the diamond (white dot Fig. [2jb)) for 
slightly different gate voltages. The rate decreases expo¬ 
nentially with temperature, indicating a thermally acti¬ 
vated process, until a background rate as described before 
is reached, giving an upper limit to the electron temper¬ 
ature. The exponential decay is potentially useful for 
thermometry [28], and again suggests involvement of the 
reservoirs with their exponential tails of the Fermi-Dirac 
distributions. At elevated temperatures, the switching 
rate exceeds the sensor bandwidth, rendering the switch¬ 
ing diamond invisible, as observed in the experiment. 

We note that the metastability shifts together with the 
DD triple points upon gate voltage adjustments. In addi¬ 
tion, it disappears when introducing significant interdot 
tunneling, and is visible also at higher vertices, see red 
circles in Fig. [ljc). This confirms the DD itself as the 
source of the switching, rather than accidental charge 
traps [29H34] . Also, the switching diamond does not ex¬ 
hibit gate-hysteresis, ruling out latching effects. Further, 
the switching persists when reducing the sensor bias to 
5 /iV, far below the A ~ 150 /ieV of the short diamond 
axis. This rules out sensor back action, which is ob¬ 
served at larger biases and exhibits a pronounced bias 
dependence [351140] . unlike the metastability seen here. 
Thus, this suggests an intrinsic DD process rather than 
an extrinsic effect as the origin for the switching. 

Based on these observations, we propose a model of 
thermal electron exchange. Inside the diamond, both 
DD one-electron levels lie below the reservoir Fermi level 
£f = 0, see Fig.[4ja,b). From a (0,1) initial state (ground 
state), the electron tunnels into a thermally activated 
hole at the same energy in the adjacent reservoir, bring¬ 
ing the DD into (0, 0). This is very slow since such a state 
lies in the exponentially suppressed tail of the Fermi func¬ 
tion, setting the overall time scale for the slow switching 
rates. Then, an electron can either fill the initial state 
(0,1) and restart the process, or it can end up in (1,0). 
Both of these transitions occur quickly - with essentially 
bare tunnel rates - since the reservoir states at the cor¬ 
responding dot energies are fully occupied. Given neg¬ 
ligible interdot tunneling, the (1,0) state is metastable, 
i.e. long lived, and is detected by the charge sensor. Al¬ 
ternatively, the process can go through (1,1) instead of 
(0,0) in a similar way. 

To test this model in the experiment, we note that it 
predicts the intermediate states (0,0) and (1,1) which 
were not seen so far (see Fig.[3|e)). If these are occupied 
long enough, they are detectable with the charge sensor. 
Thus, we decrease the tunnel coupling to the reservoirs, 
increasing the time spent in the intermediate states. Fur¬ 
ther, we heat the sample to 200 mK in order to obtain 
enough switches per unit time despite decreased reser¬ 
voir coupling. With these changes, real-time sensor data 
is recorded, see Fig.[4jc). Here, starting with (0,1), the 
DD goes briefly to the intermediate state (0, 0) before 
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FIG. 4: (a) Illustration of the model: I) slow, thermally ac¬ 
tivated tunneling into the reservoir, followed by II) tunneling 
into the metastable (1, 0) state, (b) Both /ii(l, 0) and ji 2 ( 0 ,1) 
but not ^ 1 , 2 ( 1 ,1) are below the reservoir Fermi energy within 
the shaded diamond of metastability, (c) Real-time sensor 
data for reduced reservoir coupling at 200 mK in the cen¬ 
ter of the diamond [white dot, Fig.[2][b)]. Switching via the 
intermediate states (0,0) and (1,1) is seen, (d) Four state 
Markov chain of the model, (e) Histogram of data from (c) 
as a function of sensor signal, normalized to the (0,1) and 
(1,0) positions, (f) For comparison, at the (0,0) triple point 
[yellow dot, Fig.[3Ib)], the (1,1) state is not seen. 

switching to (1,0). The DD is also sometimes seen in the 
(1,1) intermediate state, e.g. after t = 7s, here returning 
to where it was before the excursion to (1,1). Thus, we 
indeed observe the intermediate states as predicted. 

In the following, we provide an extension of the ortho¬ 
dox theory for transport in DDs m incorporating charge 
fluctuations. We note that the switching rates are slow 
enough to validate a semi-classical description with defi¬ 
nite occupation numbers Nj of dots j = 1,2. The charge 
fluctuations lead to switching between different configu¬ 
rations x = (ATi, 7V 2 ), which we express through a master 
equation for the occupation probabilities P x , 

dtPx = \.Px'^x'^x ~ P x T x ->x'] , ( 1 ) 

x'^x 

with Y x '^x the tunneling rate from configuration x' to 
x. In accordance with the experiment, we maintain only 
x = (1, 0), (0,1), (0, 0), (1,1) and neglect direct inter¬ 
dot tunneling. Hence, we keep only the rates (1,0) 
{(0,0), (1,1)} (0,1) for Eq. ([!]), see the Markov chain 

in Fig. B d )- By the Pauli principle, the bare tunneling 
rate Tj between dot j = 1,2 and its neighboring lead 
is weighted by the number of occupied lead states when 
tunneling onto the dot, /(/ij(7Vi, iV 2 )), and by the mim- 
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FIG. 5: Modeled quantities analogous to Fig. §a-d) (T e = 
60mK, Ti 5 2 = 20kHz), agreeing well with experiments. Tl,r 
are shown with color scale saturated at the lowest rates seen. 

ber of unoccupied lead states when tunneling out of the 
dot, 1 — N 2 )). Here /Xj(JVi, AT 2 ) is the chemical 

potential of dot j [27] . f(s) = [1 + exp (s/ksT)] -1 the 
Fermi function (with Boltzmann constant &#), and we 
have chosen the zero of energy at the Fermi level Sf = 0 
of the unbiased leads. The energy dependence of Ti^ 
mum is very weak compared to the energy dependence 
of the Fermi functions relevant here. Assuming energy 
independent T i ? 2 , this leads to the following set of rates: 


r ( o,o)- ( o,i) = r 2 /( M2 (o,i)) (2) 

r ( o,i)-Ko,o) = r 2 [i -/(M 2 (o,i))] (3) 

r ( o,i)-Ki,i) = r i/(Mi(l,l)) (4) 

r ( i,i)-K0,D = rqi - /(^(l,!))] (5) 


The remaining rates are obtained by (0,1) —► (1,0) 
and exchanging indices 1 f> 2. The stationary solution 
dtPx = 0 of Eq. Q becomes the inversion of a 4 x 4 ma¬ 
trix and leads to the results shown in Fig. [5] (a-d), which 
reproduce the main features in Fig. [3] The Fermi func¬ 
tions give the exponentially increasing tunnel rates when 
approaching the diamond boundaries. The zero detuning 
line exhibits no special feature in as also observed, 

since direct interdot tunneling is absent. All these fea¬ 
tures reproduce the main experimental observations. 

The background rates are not captured by the model, 
since it does not include photon absorption, interdot tun¬ 
neling and cotunneling. Interestingly, Pr ^0.5 exhibits 
a slight S-shape in the experiment, see Fig.|3^d), while 
the model calculates a straight line, see Fig^d). We 
note that the S-shape approaches a straight line upon 
reducing sensor bias. Thus, this is a sensor back action 
effect m not included in the model. We emphasize that 
only this S-shape is a back action effect - the exchange 


process itself is present in absence of charge sensing. 

As a final test, we perform a quantitative analysis of 
real-time data as shown in Fig.[4^c). We prepare his¬ 
tograms displaying the total time spent in each state. 
Close to the lower triple point, the average time spent 
in (0, 0) becomes of the same order as the time spent in 
(0,1) and (1,0), while (1,1) is almost never populated 
[44] . as seen in Fig.[4^f). In the center of the diamond, 
the DD spends most of its time in (0,1) and (1,0), and 
is rarely in (0,0) or (1,1), see Fig.^e). Pairs of simi¬ 
lar height peaks result when T i ^ T 2 . The ratio of the 
large to small peak height is given by a Boltzmann factor 
exp(A^/(fcsTe)), where Ae is the energy difference be¬ 
tween DD levels and reservoir. With T e = 200 mK, good 
agreement is found with Ae obtained using the lever arm 
extracted from fitting a Fermi function to the reservoir 
transitions. Thus, this again confirms the model. 

In summary, we report intrinsic metastable charge 
state switching within diamond shaped regions in a DD. 
The thermally activated fluctuations involve a fast elec¬ 
tron exchange with the leads, leading to an apparent tun¬ 
neling between left and right dot when direct interdot 
tunneling is negligible. An extended theory explains the 
observations very well and predicts intermediate charge 
states which are observed when reducing the reservoir 
tunnel rates below the sensor bandwidth. We emphasize 
that such thermally activated exchange continues to oc¬ 
cur when interdot tunneling is present - i.e. in absence 
of metastability - or even outside the diamond region. 
In these cases, it is directly detectable via intermediate 
states only when the sensor bandwidth exceeds the bare 
reservoir tunnel rates, but is invisible otherwise, though 
it continues to occur, possibly limiting T\. 

The exchange of a dot electron with a reservoir electron 
randomizes the DD state, thus setting an upper limit to 
Ti for both charge and spin qubits. Note that such elec¬ 
tron reservoir exchange appears also at other vertices, e.g. 
(1,1) — (0, 2) with singlet-triplet states, persists in a mag¬ 
netic field, and is generically present in any single, dou¬ 
ble, or multiple dot coupled to a reservoir, irrespective 
of the host material. Taking typical values A 5 = A/2 = 
75 /ieV, r h2 = 20 MHz, T e = 100 mK, one obtains an up¬ 
per bound Ti = T -1 exp(Ae/^T)) ^ 0.3 ms at the cen¬ 
ter of the diamond. Shifting the occupied state towards 
the reservoir Fermi energy exponentially enhances the ex¬ 
change process, facilitating fast gate-controlled spin ini¬ 
tialization [71. Decreasing T linearly extends while 
decreasing temperature or increasing the energy splitting 
As using gate voltages does so exponentially, until the 
spin-phonon coupling dominates or another process such 
as the background rate or sensor back action limits T\. 
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